

%% maps of districts, regions, PRI-PVEM coalition configuration, and campaign spending

% load shapefiles
load('../raw_data/shapefiles.mat')


%% map of districts and regions

hold on
% Region 1
mapspec = makesymbolspec('Polygon',{'Default','FaceColor',[0.2,0.2,0.2]});
mapshow(BC,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Chih,'SymbolSpec',mapspec)
mapshow(Dgo,'SymbolSpec',mapspec)
mapshow(Jal,'SymbolSpec',mapspec)
mapshow(Nay,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Sin,'SymbolSpec',mapspec)
mapshow(Son,'SymbolSpec',mapspec)
% Region 2
mapspec = makesymbolspec('Polygon',{'Default','FaceColor',[0.4,0.4,0.4]});
mapshow(Ags,'SymbolSpec',mapspec)
mapshow(Coah,'SymbolSpec',mapspec)
mapshow(Gto,'SymbolSpec',mapspec)
mapshow(NL,'SymbolSpec',mapspec)
mapshow(Qro,'SymbolSpec',mapspec)
mapshow(SLP,'SymbolSpec',mapspec)
mapshow(Tamps,'SymbolSpec',mapspec)
mapshow(Zac,'SymbolSpec',mapspec)
% Region 3
mapspec = makesymbolspec('Polygon',{'Default','FaceColor',[1,1,1]});
mapshow(Camp,'SymbolSpec',mapspec)
mapshow(Chis,'SymbolSpec',mapspec)
mapshow(Oax,'SymbolSpec',mapspec)
mapshow(QR,'SymbolSpec',mapspec)
mapshow(Tab,'SymbolSpec',mapspec)
mapshow(Ver,'SymbolSpec',mapspec)
mapshow(Yuc,'SymbolSpec',mapspec)
% Region 4
mapspec = makesymbolspec('Polygon',{'Default','FaceColor',[0.6,0.6,0.6]});
mapshow(DF,'SymbolSpec',mapspec)
mapshow(Gro,'SymbolSpec',mapspec)
mapshow(Mor,'SymbolSpec',mapspec)
mapshow(Pue,'SymbolSpec',mapspec)
mapshow(Tlax,'SymbolSpec',mapspec)
% Region 5
mapspec = makesymbolspec('Polygon',{'Default','FaceColor',[0.8,0.8,0.8]});
mapshow(Col,'SymbolSpec',mapspec)
mapshow(Hgo,'SymbolSpec',mapspec)
mapshow(Mex,'SymbolSpec',mapspec)
mapshow(Mich,'SymbolSpec',mapspec)
axis image off manual
hold off
saveas(gcf,'../figures/figA1.png')


%% map of 2012 PRI-PVEM coalition

e12 = readtable('../raw_data/election_results_2012.csv');
for d=1:3
    Ags(d).candCM = e12.candCM(d); %#ok<*SAGROW>
end
for d=1:8
    BC(d).candCM = e12.candCM(3+d);
end
for d=1:2
    BCS(d).candCM = e12.candCM(11+d);
end
for d=1:2
    Camp(d).candCM = e12.candCM(13+d);
end
for d=1:7
    Coah(d).candCM = e12.candCM(15+d);
end
for d=1:2
    Col(d).candCM = e12.candCM(22+d);
end
for d=1:12
    Chis(d).candCM = e12.candCM(24+d);
end
for d=1:9
    Chih(d).candCM = e12.candCM(36+d);
end
for d=1:27
    DF(d).candCM = e12.candCM(45+d);
end
for d=1:4
    Dgo(d).candCM = e12.candCM(72+d);
end
for d=1:14
    Gto(d).candCM = e12.candCM(76+d);
end
for d=1:9
    Gro(d).candCM = e12.candCM(90+d);
end
for d=1:7
    Hgo(d).candCM = e12.candCM(99+d);
end
for d=1:19
    Jal(d).candCM = e12.candCM(106+d);
end
for d=1:40
    Mex(d).candCM = e12.candCM(125+d);
end
for d=1:12
    Mich(d).candCM = e12.candCM(165+d);
end
for d=1:5
    Mor(d).candCM = e12.candCM(177+d);
end
for d=1:3
    Nay(d).candCM = e12.candCM(182+d);
end
for d=1:12
    NL(d).candCM = e12.candCM(185+d);
end
for d=1:11
    Oax(d).candCM = e12.candCM(197+d);
end
for d=1:16
    Pue(d).candCM = e12.candCM(208+d);
end
for d=1:4
    Qro(d).candCM = e12.candCM(224+d);
end
for d=1:3
    QR(d).candCM = e12.candCM(228+d);
end
for d=1:7
    SLP(d).candCM = e12.candCM(231+d);
end
for d=1:8
    Sin(d).candCM = e12.candCM(238+d);
end
for d=1:7
    Son(d).candCM = e12.candCM(246+d);
end
for d=1:6
    Tab(d).candCM = e12.candCM(253+d);
end
for d=1:8
    Tamps(d).candCM = e12.candCM(259+d);
end
for d=1:3
    Tlax(d).candCM = e12.candCM(267+d);
end
for d=1:21
    Ver(d).candCM = e12.candCM(270+d);
end
for d=1:5
    Yuc(d).candCM = e12.candCM(291+d);
end
for d=1:4
    Zac(d).candCM = e12.candCM(296+d);
end

hold on
% Region 1
mapspec = makesymbolspec('Polygon',{'candCM',0,'FaceColor',[1,1,1]},...
         {'candCM',1,'FaceColor',[0,0.5,0]},{'candCM',2,'FaceColor',[0,1,0]});
mapshow(BC,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Chih,'SymbolSpec',mapspec)
mapshow(Dgo,'SymbolSpec',mapspec)
mapshow(Jal,'SymbolSpec',mapspec)
mapshow(Nay,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Sin,'SymbolSpec',mapspec)
mapshow(Son,'SymbolSpec',mapspec)
% Region 2
mapshow(Ags,'SymbolSpec',mapspec)
mapshow(Coah,'SymbolSpec',mapspec)
mapshow(Gto,'SymbolSpec',mapspec)
mapshow(NL,'SymbolSpec',mapspec)
mapshow(Qro,'SymbolSpec',mapspec)
mapshow(SLP,'SymbolSpec',mapspec)
mapshow(Tamps,'SymbolSpec',mapspec)
mapshow(Zac,'SymbolSpec',mapspec)
% Region 3
mapshow(Camp,'SymbolSpec',mapspec)
mapshow(Chis,'SymbolSpec',mapspec)
mapshow(Oax,'SymbolSpec',mapspec)
mapshow(QR,'SymbolSpec',mapspec)
mapshow(Tab,'SymbolSpec',mapspec)
mapshow(Ver,'SymbolSpec',mapspec)
mapshow(Yuc,'SymbolSpec',mapspec)
% Region 4
mapshow(DF,'SymbolSpec',mapspec)
mapshow(Gro,'SymbolSpec',mapspec)
mapshow(Mor,'SymbolSpec',mapspec)
mapshow(Pue,'SymbolSpec',mapspec)
mapshow(Tlax,'SymbolSpec',mapspec)
% Region 5
mapshow(Col,'SymbolSpec',mapspec)
mapshow(Hgo,'SymbolSpec',mapspec)
mapshow(Mex,'SymbolSpec',mapspec)
mapshow(Mich,'SymbolSpec',mapspec)
axis image off manual
hold off
saveas(gcf,'../figures/fig1.png')


%% campaign spending maps

cs12 = readtable('../raw_data/campaign_spending_2012.csv');

% MP
for d=1:3
    Ags(d).CS = cs12.MP(d);
end
for d=1:8
    BC(d).CS = cs12.MP(3+d);
end
for d=1:2
    BCS(d).CS = cs12.MP(11+d);
end
for d=1:2
    Camp(d).CS = cs12.MP(13+d);
end
for d=1:7
    Coah(d).CS = cs12.MP(15+d);
end
for d=1:2
    Col(d).CS = cs12.MP(22+d);
end
for d=1:12
    Chis(d).CS = cs12.MP(24+d);
end
for d=1:9
    Chih(d).CS = cs12.MP(36+d);
end
for d=1:27
    DF(d).CS = cs12.MP(45+d);
end
for d=1:4
    Dgo(d).CS = cs12.MP(72+d);
end
for d=1:14
    Gto(d).CS = cs12.MP(76+d);
end
for d=1:9
    Gro(d).CS = cs12.MP(90+d);
end
for d=1:7
    Hgo(d).CS = cs12.MP(99+d);
end
for d=1:19
    Jal(d).CS = cs12.MP(106+d);
end
for d=1:40
    Mex(d).CS = cs12.MP(125+d);
end
for d=1:12
    Mich(d).CS = cs12.MP(165+d);
end
for d=1:5
    Mor(d).CS = cs12.MP(177+d);
end
for d=1:3
    Nay(d).CS = cs12.MP(182+d);
end
for d=1:12
    NL(d).CS = cs12.MP(185+d);
end
for d=1:11
    Oax(d).CS = cs12.MP(197+d);
end
for d=1:16
    Pue(d).CS = cs12.MP(208+d);
end
for d=1:4
    Qro(d).CS = cs12.MP(224+d);
end
for d=1:3
    QR(d).CS = cs12.MP(228+d);
end
for d=1:7
    SLP(d).CS = cs12.MP(231+d);
end
for d=1:8
    Sin(d).CS = cs12.MP(238+d);
end
for d=1:7
    Son(d).CS = cs12.MP(246+d);
end
for d=1:6
    Tab(d).CS = cs12.MP(253+d);
end
for d=1:8
    Tamps(d).CS = cs12.MP(259+d);
end
for d=1:3
    Tlax(d).CS = cs12.MP(267+d);
end
for d=1:21
    Ver(d).CS = cs12.MP(270+d);
end
for d=1:5
    Yuc(d).CS = cs12.MP(291+d);
end
for d=1:4
    Zac(d).CS = cs12.MP(296+d);
end

hold on
% Region 1
mapspec = makesymbolspec('Polygon',{'CS',[prctile(cs12.MP,0) prctile(cs12.MP,20)],'FaceColor',[0.2,0.2,0.2]},...
    {'CS',[prctile(cs12.MP,20)+10^(-10) prctile(cs12.MP,40)],'FaceColor',[0.35,0.35,0.35]},...
    {'CS',[prctile(cs12.MP,40)+10^(-10) prctile(cs12.MP,60)],'FaceColor',[0.5,0.5,0.5]},...
    {'CS',[prctile(cs12.MP,60)+10^(-10) prctile(cs12.MP,80)],'FaceColor',[0.7,0.7,0.7]},...
    {'CS',[prctile(cs12.MP,80)+10^(-10) prctile(cs12.MP,100)],'FaceColor',[0.8,0.8,0.8]});
mapshow(BC,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Chih,'SymbolSpec',mapspec)
mapshow(Dgo,'SymbolSpec',mapspec)
mapshow(Jal,'SymbolSpec',mapspec)
mapshow(Nay,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Sin,'SymbolSpec',mapspec)
mapshow(Son,'SymbolSpec',mapspec)
% Region 2
mapshow(Ags,'SymbolSpec',mapspec)
mapshow(Coah,'SymbolSpec',mapspec)
mapshow(Gto,'SymbolSpec',mapspec)
mapshow(NL,'SymbolSpec',mapspec)
mapshow(Qro,'SymbolSpec',mapspec)
mapshow(SLP,'SymbolSpec',mapspec)
mapshow(Tamps,'SymbolSpec',mapspec)
mapshow(Zac,'SymbolSpec',mapspec)
% Region 3
mapshow(Camp,'SymbolSpec',mapspec)
mapshow(Chis,'SymbolSpec',mapspec)
mapshow(Oax,'SymbolSpec',mapspec)
mapshow(QR,'SymbolSpec',mapspec)
mapshow(Tab,'SymbolSpec',mapspec)
mapshow(Ver,'SymbolSpec',mapspec)
mapshow(Yuc,'SymbolSpec',mapspec)
% Region 4
mapshow(DF,'SymbolSpec',mapspec)
mapshow(Gro,'SymbolSpec',mapspec)
mapshow(Mor,'SymbolSpec',mapspec)
mapshow(Pue,'SymbolSpec',mapspec)
mapshow(Tlax,'SymbolSpec',mapspec)
% Region 5
mapshow(Col,'SymbolSpec',mapspec)
mapshow(Hgo,'SymbolSpec',mapspec)
mapshow(Mex,'SymbolSpec',mapspec)
mapshow(Mich,'SymbolSpec',mapspec)
axis image off manual
hold off
saveas(gcf,'../figures/figA3a.png')

% NA
for d=1:3
    Ags(d).CS = cs12.NA(d);
end
for d=1:8
    BC(d).CS = cs12.NA(3+d);
end
for d=1:2
    BCS(d).CS = cs12.NA(11+d);
end
for d=1:2
    Camp(d).CS = cs12.NA(13+d);
end
for d=1:7
    Coah(d).CS = cs12.NA(15+d);
end
for d=1:2
    Col(d).CS = cs12.NA(22+d);
end
for d=1:12
    Chis(d).CS = cs12.NA(24+d);
end
for d=1:9
    Chih(d).CS = cs12.NA(36+d);
end
for d=1:27
    DF(d).CS = cs12.NA(45+d);
end
for d=1:4
    Dgo(d).CS = cs12.NA(72+d);
end
for d=1:14
    Gto(d).CS = cs12.NA(76+d);
end
for d=1:9
    Gro(d).CS = cs12.NA(90+d);
end
for d=1:7
    Hgo(d).CS = cs12.NA(99+d);
end
for d=1:19
    Jal(d).CS = cs12.NA(106+d);
end
for d=1:40
    Mex(d).CS = cs12.NA(125+d);
end
for d=1:12
    Mich(d).CS = cs12.NA(165+d);
end
for d=1:5
    Mor(d).CS = cs12.NA(177+d);
end
for d=1:3
    Nay(d).CS = cs12.NA(182+d);
end
for d=1:12
    NL(d).CS = cs12.NA(185+d);
end
for d=1:11
    Oax(d).CS = cs12.NA(197+d);
end
for d=1:16
    Pue(d).CS = cs12.NA(208+d);
end
for d=1:4
    Qro(d).CS = cs12.NA(224+d);
end
for d=1:3
    QR(d).CS = cs12.NA(228+d);
end
for d=1:7
    SLP(d).CS = cs12.NA(231+d);
end
for d=1:8
    Sin(d).CS = cs12.NA(238+d);
end
for d=1:7
    Son(d).CS = cs12.NA(246+d);
end
for d=1:6
    Tab(d).CS = cs12.NA(253+d);
end
for d=1:8
    Tamps(d).CS = cs12.NA(259+d);
end
for d=1:3
    Tlax(d).CS = cs12.NA(267+d);
end
for d=1:21
    Ver(d).CS = cs12.NA(270+d);
end
for d=1:5
    Yuc(d).CS = cs12.NA(291+d);
end
for d=1:4
    Zac(d).CS = cs12.NA(296+d);
end

hold on
% Region 1
mapspec = makesymbolspec('Polygon',{'CS',[prctile(cs12.NA,0) prctile(cs12.NA,20)],'FaceColor',[0.2,0.2,0.2]},...
    {'CS',[prctile(cs12.NA,20)+10^(-10) prctile(cs12.NA,40)],'FaceColor',[0.35,0.35,0.35]},...
    {'CS',[prctile(cs12.NA,40)+10^(-10) prctile(cs12.NA,60)],'FaceColor',[0.5,0.5,0.5]},...
    {'CS',[prctile(cs12.NA,60)+10^(-10) prctile(cs12.NA,80)],'FaceColor',[0.7,0.7,0.7]},...
    {'CS',[prctile(cs12.NA,80)+10^(-10) prctile(cs12.NA,100)],'FaceColor',[0.8,0.8,0.8]});
mapshow(BC,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Chih,'SymbolSpec',mapspec)
mapshow(Dgo,'SymbolSpec',mapspec)
mapshow(Jal,'SymbolSpec',mapspec)
mapshow(Nay,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Sin,'SymbolSpec',mapspec)
mapshow(Son,'SymbolSpec',mapspec)
% Region 2
mapshow(Ags,'SymbolSpec',mapspec)
mapshow(Coah,'SymbolSpec',mapspec)
mapshow(Gto,'SymbolSpec',mapspec)
mapshow(NL,'SymbolSpec',mapspec)
mapshow(Qro,'SymbolSpec',mapspec)
mapshow(SLP,'SymbolSpec',mapspec)
mapshow(Tamps,'SymbolSpec',mapspec)
mapshow(Zac,'SymbolSpec',mapspec)
% Region 3
mapshow(Camp,'SymbolSpec',mapspec)
mapshow(Chis,'SymbolSpec',mapspec)
mapshow(Oax,'SymbolSpec',mapspec)
mapshow(QR,'SymbolSpec',mapspec)
mapshow(Tab,'SymbolSpec',mapspec)
mapshow(Ver,'SymbolSpec',mapspec)
mapshow(Yuc,'SymbolSpec',mapspec)
% Region 4
mapshow(DF,'SymbolSpec',mapspec)
mapshow(Gro,'SymbolSpec',mapspec)
mapshow(Mor,'SymbolSpec',mapspec)
mapshow(Pue,'SymbolSpec',mapspec)
mapshow(Tlax,'SymbolSpec',mapspec)
% Region 5
mapshow(Col,'SymbolSpec',mapspec)
mapshow(Hgo,'SymbolSpec',mapspec)
mapshow(Mex,'SymbolSpec',mapspec)
mapshow(Mich,'SymbolSpec',mapspec)
axis image off manual
hold off
saveas(gcf,'../figures/figA3b.png')

% PVEM (alone)
for d=1:3
    Ags(d).CS = cs12.PVEM(d);
end
for d=1:8
    BC(d).CS = cs12.PVEM(3+d);
end
for d=1:2
    BCS(d).CS = cs12.PVEM(11+d);
end
for d=1:2
    Camp(d).CS = cs12.PVEM(13+d);
end
for d=1:7
    Coah(d).CS = cs12.PVEM(15+d);
end
for d=1:2
    Col(d).CS = cs12.PVEM(22+d);
end
for d=1:12
    Chis(d).CS = cs12.PVEM(24+d);
end
for d=1:9
    Chih(d).CS = cs12.PVEM(36+d);
end
for d=1:27
    DF(d).CS = cs12.PVEM(45+d);
end
for d=1:4
    Dgo(d).CS = cs12.PVEM(72+d);
end
for d=1:14
    Gto(d).CS = cs12.PVEM(76+d);
end
for d=1:9
    Gro(d).CS = cs12.PVEM(90+d);
end
for d=1:7
    Hgo(d).CS = cs12.PVEM(99+d);
end
for d=1:19
    Jal(d).CS = cs12.PVEM(106+d);
end
for d=1:40
    Mex(d).CS = cs12.PVEM(125+d);
end
for d=1:12
    Mich(d).CS = cs12.PVEM(165+d);
end
for d=1:5
    Mor(d).CS = cs12.PVEM(177+d);
end
for d=1:3
    Nay(d).CS = cs12.PVEM(182+d);
end
for d=1:12
    NL(d).CS = cs12.PVEM(185+d);
end
for d=1:11
    Oax(d).CS = cs12.PVEM(197+d);
end
for d=1:16
    Pue(d).CS = cs12.PVEM(208+d);
end
for d=1:4
    Qro(d).CS = cs12.PVEM(224+d);
end
for d=1:3
    QR(d).CS = cs12.PVEM(228+d);
end
for d=1:7
    SLP(d).CS = cs12.PVEM(231+d);
end
for d=1:8
    Sin(d).CS = cs12.PVEM(238+d);
end
for d=1:7
    Son(d).CS = cs12.PVEM(246+d);
end
for d=1:6
    Tab(d).CS = cs12.PVEM(253+d);
end
for d=1:8
    Tamps(d).CS = cs12.PVEM(259+d);
end
for d=1:3
    Tlax(d).CS = cs12.PVEM(267+d);
end
for d=1:21
    Ver(d).CS = cs12.PVEM(270+d);
end
for d=1:5
    Yuc(d).CS = cs12.PVEM(291+d);
end
for d=1:4
    Zac(d).CS = cs12.PVEM(296+d);
end

hold on
% Region 1
mapspec = makesymbolspec('Polygon',{'CS',[prctile(cs12.PVEM(cs12.PVEM>0),0) prctile(cs12.PVEM(cs12.PVEM>0),20)],'FaceColor',[0.2,0.2,0.2]},...
    {'CS',[prctile(cs12.PVEM(cs12.PVEM>0),20)+10^(-10) prctile(cs12.PVEM(cs12.PVEM>0),40)],'FaceColor',[0.35,0.35,0.35]},...
    {'CS',[prctile(cs12.PVEM(cs12.PVEM>0),40)+10^(-10) prctile(cs12.PVEM(cs12.PVEM>0),60)],'FaceColor',[0.5,0.5,0.5]},...
    {'CS',[prctile(cs12.PVEM(cs12.PVEM>0),60)+10^(-10) prctile(cs12.PVEM(cs12.PVEM>0),80)],'FaceColor',[0.7,0.7,0.7]},...
    {'CS',[prctile(cs12.PVEM(cs12.PVEM>0),80)+10^(-10) prctile(cs12.PVEM(cs12.PVEM>0),100)],'FaceColor',[0.8,0.8,0.8]},...
    {'CS',0,'FaceColor',[1,1,1]});
mapshow(BC,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Chih,'SymbolSpec',mapspec)
mapshow(Dgo,'SymbolSpec',mapspec)
mapshow(Jal,'SymbolSpec',mapspec)
mapshow(Nay,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Sin,'SymbolSpec',mapspec)
mapshow(Son,'SymbolSpec',mapspec)
% Region 2
mapshow(Ags,'SymbolSpec',mapspec)
mapshow(Coah,'SymbolSpec',mapspec)
mapshow(Gto,'SymbolSpec',mapspec)
mapshow(NL,'SymbolSpec',mapspec)
mapshow(Qro,'SymbolSpec',mapspec)
mapshow(SLP,'SymbolSpec',mapspec)
mapshow(Tamps,'SymbolSpec',mapspec)
mapshow(Zac,'SymbolSpec',mapspec)
% Region 3
mapshow(Camp,'SymbolSpec',mapspec)
mapshow(Chis,'SymbolSpec',mapspec)
mapshow(Oax,'SymbolSpec',mapspec)
mapshow(QR,'SymbolSpec',mapspec)
mapshow(Tab,'SymbolSpec',mapspec)
mapshow(Ver,'SymbolSpec',mapspec)
mapshow(Yuc,'SymbolSpec',mapspec)
% Region 4
mapshow(DF,'SymbolSpec',mapspec)
mapshow(Gro,'SymbolSpec',mapspec)
mapshow(Mor,'SymbolSpec',mapspec)
mapshow(Pue,'SymbolSpec',mapspec)
mapshow(Tlax,'SymbolSpec',mapspec)
% Region 5
mapshow(Col,'SymbolSpec',mapspec)
mapshow(Hgo,'SymbolSpec',mapspec)
mapshow(Mex,'SymbolSpec',mapspec)
mapshow(Mich,'SymbolSpec',mapspec)
axis image off manual
hold off
saveas(gcf,'../figures/figA3c.png')

% PRI (alone)
for d=1:3
    Ags(d).CS = cs12.PRI(d);
end
for d=1:8
    BC(d).CS = cs12.PRI(3+d);
end
for d=1:2
    BCS(d).CS = cs12.PRI(11+d);
end
for d=1:2
    Camp(d).CS = cs12.PRI(13+d);
end
for d=1:7
    Coah(d).CS = cs12.PRI(15+d);
end
for d=1:2
    Col(d).CS = cs12.PRI(22+d);
end
for d=1:12
    Chis(d).CS = cs12.PRI(24+d);
end
for d=1:9
    Chih(d).CS = cs12.PRI(36+d);
end
for d=1:27
    DF(d).CS = cs12.PRI(45+d);
end
for d=1:4
    Dgo(d).CS = cs12.PRI(72+d);
end
for d=1:14
    Gto(d).CS = cs12.PRI(76+d);
end
for d=1:9
    Gro(d).CS = cs12.PRI(90+d);
end
for d=1:7
    Hgo(d).CS = cs12.PRI(99+d);
end
for d=1:19
    Jal(d).CS = cs12.PRI(106+d);
end
for d=1:40
    Mex(d).CS = cs12.PRI(125+d);
end
for d=1:12
    Mich(d).CS = cs12.PRI(165+d);
end
for d=1:5
    Mor(d).CS = cs12.PRI(177+d);
end
for d=1:3
    Nay(d).CS = cs12.PRI(182+d);
end
for d=1:12
    NL(d).CS = cs12.PRI(185+d);
end
for d=1:11
    Oax(d).CS = cs12.PRI(197+d);
end
for d=1:16
    Pue(d).CS = cs12.PRI(208+d);
end
for d=1:4
    Qro(d).CS = cs12.PRI(224+d);
end
for d=1:3
    QR(d).CS = cs12.PRI(228+d);
end
for d=1:7
    SLP(d).CS = cs12.PRI(231+d);
end
for d=1:8
    Sin(d).CS = cs12.PRI(238+d);
end
for d=1:7
    Son(d).CS = cs12.PRI(246+d);
end
for d=1:6
    Tab(d).CS = cs12.PRI(253+d);
end
for d=1:8
    Tamps(d).CS = cs12.PRI(259+d);
end
for d=1:3
    Tlax(d).CS = cs12.PRI(267+d);
end
for d=1:21
    Ver(d).CS = cs12.PRI(270+d);
end
for d=1:5
    Yuc(d).CS = cs12.PRI(291+d);
end
for d=1:4
    Zac(d).CS = cs12.PRI(296+d);
end

hold on
% Region 1
mapspec = makesymbolspec('Polygon',{'CS',[prctile(cs12.PRI(cs12.PRI>0),0) prctile(cs12.PRI(cs12.PRI>0),20)],'FaceColor',[0.2,0.2,0.2]},...
    {'CS',[prctile(cs12.PRI(cs12.PRI>0),20)+10^(-10) prctile(cs12.PRI(cs12.PRI>0),40)],'FaceColor',[0.35,0.35,0.35]},...
    {'CS',[prctile(cs12.PRI(cs12.PRI>0),40)+10^(-10) prctile(cs12.PRI(cs12.PRI>0),60)],'FaceColor',[0.5,0.5,0.5]},...
    {'CS',[prctile(cs12.PRI(cs12.PRI>0),60)+10^(-10) prctile(cs12.PRI(cs12.PRI>0),80)],'FaceColor',[0.7,0.7,0.7]},...
    {'CS',[prctile(cs12.PRI(cs12.PRI>0),80)+10^(-10) prctile(cs12.PRI(cs12.PRI>0),100)],'FaceColor',[0.8,0.8,0.8]},...
    {'CS',0,'FaceColor',[1,1,1]});
mapshow(BC,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Chih,'SymbolSpec',mapspec)
mapshow(Dgo,'SymbolSpec',mapspec)
mapshow(Jal,'SymbolSpec',mapspec)
mapshow(Nay,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Sin,'SymbolSpec',mapspec)
mapshow(Son,'SymbolSpec',mapspec)
% Region 2
mapshow(Ags,'SymbolSpec',mapspec)
mapshow(Coah,'SymbolSpec',mapspec)
mapshow(Gto,'SymbolSpec',mapspec)
mapshow(NL,'SymbolSpec',mapspec)
mapshow(Qro,'SymbolSpec',mapspec)
mapshow(SLP,'SymbolSpec',mapspec)
mapshow(Tamps,'SymbolSpec',mapspec)
mapshow(Zac,'SymbolSpec',mapspec)
% Region 3
mapshow(Camp,'SymbolSpec',mapspec)
mapshow(Chis,'SymbolSpec',mapspec)
mapshow(Oax,'SymbolSpec',mapspec)
mapshow(QR,'SymbolSpec',mapspec)
mapshow(Tab,'SymbolSpec',mapspec)
mapshow(Ver,'SymbolSpec',mapspec)
mapshow(Yuc,'SymbolSpec',mapspec)
% Region 4
mapshow(DF,'SymbolSpec',mapspec)
mapshow(Gro,'SymbolSpec',mapspec)
mapshow(Mor,'SymbolSpec',mapspec)
mapshow(Pue,'SymbolSpec',mapspec)
mapshow(Tlax,'SymbolSpec',mapspec)
% Region 5
mapshow(Col,'SymbolSpec',mapspec)
mapshow(Hgo,'SymbolSpec',mapspec)
mapshow(Mex,'SymbolSpec',mapspec)
mapshow(Mich,'SymbolSpec',mapspec)
axis image off manual
hold off
saveas(gcf,'../figures/figA3d.png')

% CM
for d=1:3
    Ags(d).CS = cs12.CM(d);
end
for d=1:8
    BC(d).CS = cs12.CM(3+d);
end
for d=1:2
    BCS(d).CS = cs12.CM(11+d);
end
for d=1:2
    Camp(d).CS = cs12.CM(13+d);
end
for d=1:7
    Coah(d).CS = cs12.CM(15+d);
end
for d=1:2
    Col(d).CS = cs12.CM(22+d);
end
for d=1:12
    Chis(d).CS = cs12.CM(24+d);
end
for d=1:9
    Chih(d).CS = cs12.CM(36+d);
end
for d=1:27
    DF(d).CS = cs12.CM(45+d);
end
for d=1:4
    Dgo(d).CS = cs12.CM(72+d);
end
for d=1:14
    Gto(d).CS = cs12.CM(76+d);
end
for d=1:9
    Gro(d).CS = cs12.CM(90+d);
end
for d=1:7
    Hgo(d).CS = cs12.CM(99+d);
end
for d=1:19
    Jal(d).CS = cs12.CM(106+d);
end
for d=1:40
    Mex(d).CS = cs12.CM(125+d);
end
for d=1:12
    Mich(d).CS = cs12.CM(165+d);
end
for d=1:5
    Mor(d).CS = cs12.CM(177+d);
end
for d=1:3
    Nay(d).CS = cs12.CM(182+d);
end
for d=1:12
    NL(d).CS = cs12.CM(185+d);
end
for d=1:11
    Oax(d).CS = cs12.CM(197+d);
end
for d=1:16
    Pue(d).CS = cs12.CM(208+d);
end
for d=1:4
    Qro(d).CS = cs12.CM(224+d);
end
for d=1:3
    QR(d).CS = cs12.CM(228+d);
end
for d=1:7
    SLP(d).CS = cs12.CM(231+d);
end
for d=1:8
    Sin(d).CS = cs12.CM(238+d);
end
for d=1:7
    Son(d).CS = cs12.CM(246+d);
end
for d=1:6
    Tab(d).CS = cs12.CM(253+d);
end
for d=1:8
    Tamps(d).CS = cs12.CM(259+d);
end
for d=1:3
    Tlax(d).CS = cs12.CM(267+d);
end
for d=1:21
    Ver(d).CS = cs12.CM(270+d);
end
for d=1:5
    Yuc(d).CS = cs12.CM(291+d);
end
for d=1:4
    Zac(d).CS = cs12.CM(296+d);
end

hold on
% Region 1
mapspec = makesymbolspec('Polygon',{'CS',[prctile(cs12.CM(cs12.CM>0),0) prctile(cs12.CM(cs12.CM>0),20)],'FaceColor',[0.2,0.2,0.2]},...
    {'CS',[prctile(cs12.CM(cs12.CM>0),20)+10^(-10) prctile(cs12.CM(cs12.CM>0),40)],'FaceColor',[0.35,0.35,0.35]},...
    {'CS',[prctile(cs12.CM(cs12.CM>0),40)+10^(-10) prctile(cs12.CM(cs12.CM>0),60)],'FaceColor',[0.5,0.5,0.5]},...
    {'CS',[prctile(cs12.CM(cs12.CM>0),60)+10^(-10) prctile(cs12.CM(cs12.CM>0),80)],'FaceColor',[0.7,0.7,0.7]},...
    {'CS',[prctile(cs12.CM(cs12.CM>0),80)+10^(-10) prctile(cs12.CM(cs12.CM>0),100)],'FaceColor',[0.8,0.8,0.8]},...
    {'CS',0,'FaceColor',[1,1,1]});
mapshow(BC,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Chih,'SymbolSpec',mapspec)
mapshow(Dgo,'SymbolSpec',mapspec)
mapshow(Jal,'SymbolSpec',mapspec)
mapshow(Nay,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Sin,'SymbolSpec',mapspec)
mapshow(Son,'SymbolSpec',mapspec)
% Region 2
mapshow(Ags,'SymbolSpec',mapspec)
mapshow(Coah,'SymbolSpec',mapspec)
mapshow(Gto,'SymbolSpec',mapspec)
mapshow(NL,'SymbolSpec',mapspec)
mapshow(Qro,'SymbolSpec',mapspec)
mapshow(SLP,'SymbolSpec',mapspec)
mapshow(Tamps,'SymbolSpec',mapspec)
mapshow(Zac,'SymbolSpec',mapspec)
% Region 3
mapshow(Camp,'SymbolSpec',mapspec)
mapshow(Chis,'SymbolSpec',mapspec)
mapshow(Oax,'SymbolSpec',mapspec)
mapshow(QR,'SymbolSpec',mapspec)
mapshow(Tab,'SymbolSpec',mapspec)
mapshow(Ver,'SymbolSpec',mapspec)
mapshow(Yuc,'SymbolSpec',mapspec)
% Region 4
mapshow(DF,'SymbolSpec',mapspec)
mapshow(Gro,'SymbolSpec',mapspec)
mapshow(Mor,'SymbolSpec',mapspec)
mapshow(Pue,'SymbolSpec',mapspec)
mapshow(Tlax,'SymbolSpec',mapspec)
% Region 5
mapshow(Col,'SymbolSpec',mapspec)
mapshow(Hgo,'SymbolSpec',mapspec)
mapshow(Mex,'SymbolSpec',mapspec)
mapshow(Mich,'SymbolSpec',mapspec)
axis image off manual
hold off
saveas(gcf,'../figures/figA3e.png')

% PAN
for d=1:3
    Ags(d).CS = cs12.PAN(d);
end
for d=1:8
    BC(d).CS = cs12.PAN(3+d);
end
for d=1:2
    BCS(d).CS = cs12.PAN(11+d);
end
for d=1:2
    Camp(d).CS = cs12.PAN(13+d);
end
for d=1:7
    Coah(d).CS = cs12.PAN(15+d);
end
for d=1:2
    Col(d).CS = cs12.PAN(22+d);
end
for d=1:12
    Chis(d).CS = cs12.PAN(24+d);
end
for d=1:9
    Chih(d).CS = cs12.PAN(36+d);
end
for d=1:27
    DF(d).CS = cs12.PAN(45+d);
end
for d=1:4
    Dgo(d).CS = cs12.PAN(72+d);
end
for d=1:14
    Gto(d).CS = cs12.PAN(76+d);
end
for d=1:9
    Gro(d).CS = cs12.PAN(90+d);
end
for d=1:7
    Hgo(d).CS = cs12.PAN(99+d);
end
for d=1:19
    Jal(d).CS = cs12.PAN(106+d);
end
for d=1:40
    Mex(d).CS = cs12.PAN(125+d);
end
for d=1:12
    Mich(d).CS = cs12.PAN(165+d);
end
for d=1:5
    Mor(d).CS = cs12.PAN(177+d);
end
for d=1:3
    Nay(d).CS = cs12.PAN(182+d);
end
for d=1:12
    NL(d).CS = cs12.PAN(185+d);
end
for d=1:11
    Oax(d).CS = cs12.PAN(197+d);
end
for d=1:16
    Pue(d).CS = cs12.PAN(208+d);
end
for d=1:4
    Qro(d).CS = cs12.PAN(224+d);
end
for d=1:3
    QR(d).CS = cs12.PAN(228+d);
end
for d=1:7
    SLP(d).CS = cs12.PAN(231+d);
end
for d=1:8
    Sin(d).CS = cs12.PAN(238+d);
end
for d=1:7
    Son(d).CS = cs12.PAN(246+d);
end
for d=1:6
    Tab(d).CS = cs12.PAN(253+d);
end
for d=1:8
    Tamps(d).CS = cs12.PAN(259+d);
end
for d=1:3
    Tlax(d).CS = cs12.PAN(267+d);
end
for d=1:21
    Ver(d).CS = cs12.PAN(270+d);
end
for d=1:5
    Yuc(d).CS = cs12.PAN(291+d);
end
for d=1:4
    Zac(d).CS = cs12.PAN(296+d);
end

hold on
% Region 1
mapspec = makesymbolspec('Polygon',{'CS',[prctile(cs12.PAN,0) prctile(cs12.PAN,20)],'FaceColor',[0.2,0.2,0.2]},...
    {'CS',[prctile(cs12.PAN,20)+10^(-10) prctile(cs12.PAN,40)],'FaceColor',[0.35,0.35,0.35]},...
    {'CS',[prctile(cs12.PAN,40)+10^(-10) prctile(cs12.PAN,60)],'FaceColor',[0.5,0.5,0.5]},...
    {'CS',[prctile(cs12.PAN,60)+10^(-10) prctile(cs12.PAN,80)],'FaceColor',[0.7,0.7,0.7]},...
    {'CS',[prctile(cs12.PAN,80)+10^(-10) prctile(cs12.PAN,100)],'FaceColor',[0.8,0.8,0.8]});
mapshow(BC,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Chih,'SymbolSpec',mapspec)
mapshow(Dgo,'SymbolSpec',mapspec)
mapshow(Jal,'SymbolSpec',mapspec)
mapshow(Nay,'SymbolSpec',mapspec)
mapshow(BCS,'SymbolSpec',mapspec)
mapshow(Sin,'SymbolSpec',mapspec)
mapshow(Son,'SymbolSpec',mapspec)
% Region 2
mapshow(Ags,'SymbolSpec',mapspec)
mapshow(Coah,'SymbolSpec',mapspec)
mapshow(Gto,'SymbolSpec',mapspec)
mapshow(NL,'SymbolSpec',mapspec)
mapshow(Qro,'SymbolSpec',mapspec)
mapshow(SLP,'SymbolSpec',mapspec)
mapshow(Tamps,'SymbolSpec',mapspec)
mapshow(Zac,'SymbolSpec',mapspec)
% Region 3
mapshow(Camp,'SymbolSpec',mapspec)
mapshow(Chis,'SymbolSpec',mapspec)
mapshow(Oax,'SymbolSpec',mapspec)
mapshow(QR,'SymbolSpec',mapspec)
mapshow(Tab,'SymbolSpec',mapspec)
mapshow(Ver,'SymbolSpec',mapspec)
mapshow(Yuc,'SymbolSpec',mapspec)
% Region 4
mapshow(DF,'SymbolSpec',mapspec)
mapshow(Gro,'SymbolSpec',mapspec)
mapshow(Mor,'SymbolSpec',mapspec)
mapshow(Pue,'SymbolSpec',mapspec)
mapshow(Tlax,'SymbolSpec',mapspec)
% Region 5
mapshow(Col,'SymbolSpec',mapspec)
mapshow(Hgo,'SymbolSpec',mapspec)
mapshow(Mex,'SymbolSpec',mapspec)
mapshow(Mich,'SymbolSpec',mapspec)
axis image off manual
hold off
saveas(gcf,'../figures/figA3f.png')



